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Planets grow in rotating disks of dust and gas around forming stars, some of which 
can subsequently collide in giant impacts after the gas component is removed from 
the disk’?. Monitoring programmes with the warm Spitzer mission have recorded 
substantial and rapid changes in mid-infrared output for several stars, interpreted 
as variations in the surface area of warm, dusty material ejected by planetary-scale 
collisions and heated by the central star: for example, NGC 2354-ID8 (refs. 4,5), HD 
166191 (ref. 6) and V488 Persei”. Here we report combined observations of the young 
(about 300 million years old), solar-like star ASASSN-21qj: an infrared brightening 
consistent with a blackbody temperature of 1,000 Kelvin and aluminosity that is 

4 percent that of the star lasting for about 1,000 days, partially overlapping in time 
with a complex and deep, wavelength-dependent optical eclipse that lasted for 
about 500 days. The optical eclipse started 2.5 years after the infrared brightening, 
implying an orbital period of at least that duration. These observations are consistent 
with a collision between two exoplanets of several to tens of Earth masses at 2-16 
astronomical units from the central star. Such an impact produces a hot, highly 
extended post-impact remnant with sufficient luminosity to explain the infrared 
observations. Transit of the impact debris, sheared by orbital motion into along 
cloud, causes the subsequent complex eclipse of the host star. 


The otherwise unexceptional star 2MASS JO8152329-3859234 under- 
went a sudden optical-dimming event in December 2021°” and was 
assigned the identifier ASASSN-21qj by the All Sky Automated Survey 
for Supernovae (ASAS-SN)'". Here we combine both optical (from 
the Las Cumbres Observatory Global Telescope Network (LCOGT)) 
and infrared (from the Wide-field Infrared Survey Explorer (WISE) 
satellite) observations of ASASSN-21qj for the years before and after 
this dimming event (Fig. 1). Optical multiband photometry shows a 
wavelength-dependent depth consistent with extinction by submi- 
crometre particles. About 900 days before the optical-dimming event, 
the ASASSN-21qj system showed a notable brightening in the infrared, 
of 0.4 magnitudes at 3.8 um (W1) and 0.8 magnitudes at 4.5 um (W2). 
Before this time, the infrared brightness was consistent with being 
purely stellar. The increased infrared fluxes indicate that, as well as 
the quiescent stellar flux, there was new emission at a temperature of 
approximately 1,000 K. Such a notable combination of observations, 


particularly the 2.5-year delay between the infrared and optical varia- 
tion, requires an explanation. 

The optical and infrared light curves (Fig. 1) provide key constraints 
onany proposed scenario. The post-brightening infrared fluxes in the 
Wland W2 WISE passbands are consistent with emission at a blackbody 
temperature of 1,000 + 100 K and this temperature is sustained, within 
error, for the remainder of our observation window, despite a declinein 
the total flux. For an emitter located at the distance of ASASSN-21qj from 
Earth, and the observed maximum luminosity of approximately 0.04 L,, 
this implies an emitting area of 0.01 AU’, equivalent to an object witha 
radius of 7 Ro, or 750 Reatn If this emission was from material—such as 
dust—passively heated by proximity to the star, then that material must 
have been generated and remained within about 0.1 AU to produce the 
observed temperature. 

One possible explanation is that we are observing two unrelated but 
coincidental phenomena: a warm-dust-generating collision within 
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Fig. 1| Optical and infrared photometry of ASASSN-21qj. a, Normalized 
optical photometry from ASAS-SN in the V-band and the g’-band. b, Fractional 
flux increase in brightness of ASASSN-21qj in both the W1and W2 bands, for 
which a value of 1.0 represents the stellar contribution alone. c, Calculated 
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NEOWISE colour temperature estimated from the photometry of the two 
bands. The colour temperature is plotted as zero when there is no infrared 
excess and is consistent with a temperature of 1,000 K while the excess is 
present. Error bars are shown at loconfidence. 
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Fig. 2| The light curve of ASASSN-21qj from several different photometric 
surveys and the derived transverse velocities. a, The eclipse depthis deeper 
for shorter wavelengths, indicating that the transiting material is dominated by 
submicrometre-sized grains. b, Transverse velocities derived from the 
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light-curve gradients. These are lower limits to the true velocity and thus imply 
that the transiting material is closer to the star than 16 AU. Error bars are shown 
at loconfidence. 
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Fig. 3 | The size and temporal evolution of a post-impact body. A post-impact 
body of a few to tens of Earth masses could be large enough to explain the 
observed increase in thermal emission from ASASSN-21qj and the subsequent 
infrared fluxes. a, The Hill radii for bodies at different semimajor axes 

(solid coloured lines) and the Bondi radii (dashed lines) assuming different 
compositions for the vapour of the post-impact body. The horizontal black 
dashed line shows the radius required to explain the observed flux at the 
inferred emission temperature of 1,000 K. b,c, The change in excess flux owing 
toa post-impact body with time for a simple model of cooling of a post-impact 


0.1 AU of the star, with a separate object transiting the star 900 days 
later. Two events that are themselves very rare occurring independently 
in one system is, however, highly improbable. A second explanation 
is that warm dust is generated close to the star and causes the optical 
transit, but this requires a fine-tuned configuration in which the star 
is optically blocked by scale-height variations in the resulting disk. 

Instead, we suggest that we are observing the aftermath of a single 
collision between super-Earths or mini-Neptunes—a so-called giant 
impact—between 2 and 16 AU from the star. These distances are deter- 
mined, respectively, by the delay between the infrared brightening and 
the optical eclipse (Fig. 1) and by gradients in the optical light curve 
(Fig. 2). In contrast to other extreme-debris-disk events in which the star 
heats the dust, we propose that the infrared emission is directly from 
the post-impact body” and that impact debris produced the optical 
transit. Giant impacts are a common occurrence in planet formation” 
and also occur during instabilities in older systems”*; this would explain 
the observations witha single event of atypethat is expected for systems 
with ages such as that of ASASSN-21qj. 

Giant impacts are one of the most energetic events that planets 
experience. For example, the kinetic energy of impacts between two 
half-Neptune-mass bodies range from 10* to 10** J, enough to vaporize 
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body for different power-law surface density profiles (b; equation (1)) and for 
different mass bodies (c). Each example had an initial radius of 7 R,, a radiative 
temperature of 1,000 K and a corresponding initial flux of 0.04 L,, as estimated 
for the observed infrared emitter. The solid lines are profiles that can have 
non-zero initial density at the initial emitting radius and the dashed lines are for 
ones that are forced to have zero density at the initial emitting radius. When not 
stated, the mass of the post-impact body was 50 Marn the power-law exponent 
is —2 and the size of the central region was 10 Reartn- 


the colliding bodies several times over. A large fraction of this energy 
is dissipated in the colliding bodies and post-impact bodies are sub- 
stantially melted and vaporized”"®, Furthermore, extreme torques 
exerted in impacts often produce rapidly rotating bodies”. Such low- 
density and rotationally flattened bodies can be hundreds of times 
larger than the pre-impact planets”, with correspondingly large radia- 
tive surfaces. 

Giant impacts produce substantial amounts of debris, typically 
around 1% of the colliding mass, that is injected into orbit around the 
host star’”°. Impact ejecta have a wide range of sizes, from submicro- 
metre dust to planetesimals of tens to hundreds of kilometres across, 
and often contain the most highly heated material” ~’. For sufficiently 
high impact velocities (>1 km s” for water ice and >8 km s” for forster- 
ite), asubstantial fraction of this material is vaporized” ”. Shearing of 
droplets and cooling and condensation of vapour produces a popula- 
tion of small dust grains and solid spherules. The size distribution of 
this fraction of the debris is uncertain, owing to difficulties in model- 
ling condensate nucleation and break-up, but previous work suggests 
that debris could range in size from submicrometre to decimetres””®. 
The wavelength-dependent eclipse suggests that the optical depth 
of the transiting dust cloud is dominated by submicrometre grains, 
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consistent with these previous estimates. Although the fading of the 
excess infrared flux occurs within 100 days of the start of the optical 
transit, we consider the timing to be coincidental because there is no 
clear correspondence between the light curves, for example, no change 
inthe (6-month cadence) infrared flux when the star dims in the optical 
wavelengths just before MJD 59500. 

The radiative flux from post-impact bodies has not been explored 
in depth. Computational-resource limitations make resolving the 
low-density outer regions and photosphere of post-impact bodies 
extremely challenging. However, preliminary simulations of impacts 
between super-Earth and mini-Neptunes place an approximate lower 
limit on the extent of post-impact bodies and show that post-impact 
bodies can extend to hundreds of Earth radii. Such an object radiating 
at about 1,000 K would produce a flux comparable with the 0.04 L, 
inferred from our observations. 

Independent of impact simulations, there are fundamental limits 
on the size of post-impact bodies from the Hill and Bondi radii. Any 
post-impact structure must lie within the Hill sphere, the distance 
within which the gravity of an object dominates over that of the star. 
Furthermore, only vapour within the Bondi radius would be bound 
to the post-impact body. Figure 3a shows the Hill radii at different 
distances from ASASSN-21qj (solid lines) and the Bondi radii for 
example gas species (dashed lines). Beyond 2.4 Au, the Hill radii of 
greater-than-Earth-mass bodies are large enough to accommodate 
a post-impact body capable of producing the required infrared flux 
(approximately 7 R,, dashed line). Heavier gases (H,O, SiO and SiO,) are 
also bound to bodies of more than a few Earth masses. A post-impact 
body ofa few Earth masses can hence theoretically produce the observed 
infrared emission, with a photosphere that is dust and/or vapour. 

A key line of evidence for direct detection of a post-impact body is the 
constant emission temperature. It is argued” that rocky post-impact 
bodies become optically thin at low pressures at which radiative loss 
drives rapid condensation of the rock vapour. The emission tempera- 
ture is then set by the liquid-vapour phase boundary and is constant 
until the post-impact body almost fully condenses”°””. The dew point 
(onset of condensation) and bubble point (onset of substantial vapori- 
zation) of material with bulk silicate Earth composition are similar 
(about 2,300 K, within roughly 100 K (refs. 20,28)) at the low pres- 
sures of the photosphere of a post-impact body. However, even small 
amounts of water (10° mole fraction) may lower the bubble and dew 
points for silicates by approximately 100 K (refs. 20,28). The emission 
temperature of bodies produced by collisions of protoplanets com- 
posed of rock and ices/volatiles could be buffered at about 1,000 K 
during early evolution. 

The temporal variation of flux froma post-impact body is controlled 
by evolution of the sizes of the bodies, which is governed by several 
competing factors”°”’, including radiative energy loss, viscous spread- 
ing and mass and angular momentum transfer by condensates. For 
silicate-dominated post-impact bodies, the high emission temperature 
means that radiative cooling dominates and the post-impact body con- 
tracts rapidly, fully condensing over years to thousands of years”°”’. 
The body we observed has a much lower emission temperature and 
contracted much more slowly. Figure 3b,c shows the evolution of con- 
ceptual post-impact bodies in the limiting case that radiative cooling 
dominates. If sufficient mass is injected into the outer regions (that 
is, less negative power-law exponents), the observed flux can remain 
constant for an initial period and then decays over the order of hun- 
dreds of days, inagreement with infrared observations. Further work 
is required to understand the structure and evolution of bodies pro- 
duced by impacts between super-Earths and mini-Neptunes of different 
compositions and identify temporal-flux variations consistent with 
our observations. 
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Methods 


The stellar properties of ASASSN-21qj (Gaia DR3 5539970601632026752 
= 2MASS J08152329-3859234) are listed in Extended Data Table 1, 
showing that ASASSN-21qj is consistent with being a G2-type dwarf 
star. Where necessary, we assume a stellar mass equal to the Sun. 
ASASSN-21qj has a neighbour (Gaia DR3 5539970597334497024 
= 2MASS J08152298-3859244), which is a visual double. On the basis 
of the Gaia DR3 mean International Celestial Reference System (ICRS) 
position for epoch 2016.0, the visual companion lies at a separation 
P =3,738.243 + 0.062 mas and at position angle 0 = 249.977°. Their par- 
allaxes (w = 1.7631 + 0.0112 mas versus 1.4711 + 0.0523 mas) differ by 5.50 
and proper motions (H, = —9.692 + 0.012, u5 = 7.349 + 0.012 mas year? 
versus H, = —0.114 + 0.055, 45 = 6.419 + 0.053 mas year”) differ by a factor 
of 2. The large differences in distance and proper motions suggest 
that these stars are not associated. 

The stellar photospheric flux was estimated by fitting stellar models 
to GAIA, APASS and DENIS and WISE optical/near-infrared photometry 
(the 2MASS/ photometry is an upper limit and H and K, are flagged 
as contaminated). Extended Data Fig. 1 shows the resulting models 
(the dashed line is discussed below). The ALLWISE photometry (around 
2010) is consistent with the Near-Earth Object Wide-field Infrared 
Survey Explorer (NEOWISE) photometry pre-brightening nearly 
10 years later. 

A fit using the method in ref. 30 with GAIA and DENIS photometry 
finds that the WISE W1/2 fluxes are about 20% too high, but better 
agreementis found with APASS photometry. The difference is explained 
by the fact that WISE and APASS have lower spatial resolution and 
include the flux of the approximately 2 mag fainter visual double to 
the west of ASASSN-21qj (whichis visible in 2MASS). The best-fit stellar 
effective temperature is 5,560 + 100 K. We donot find that reddening is 
needed for these models; although there is relatively little photometry 
with which to strongly constrain both T,,-and Ay, the Tep from GAIA DR3 
is 5,760 + 10 K, which suggests that the conclusion of little reddening is 
valid. Primarily, we conclude that there is no indication that ASASSN-21qj 
showed evidence for an infrared excess before the brightening seen 
by NEOWISE. 

To estimate the infrared excess properties, we fit the same models 
but nowincluding the first five post-brightening NEOWISE data points. 
We correct for the nearby source by using the GAIA photometry 
or the star and the post-brightening/pre-brightening difference 
for the WISE excess flux. This fit yields the fractional luminosity 
L aus/L, = 0.04 + 0.005 and a dust temperature 950 + 30 K. To estimate 
the dust temperature as a function of time, we simply subtract the 
median pre-brightening W1/2 fluxes, with the r.m.s. of these values 
as the uncertainty (see Fig. 1). The temperature uncertainty increases 
as the excess fades; the excess flux uncertainty depends on both the 
observed flux and the stellar flux, and as the excess decreases, the stellar 
uncertainty (whichis constant) becomes an increasingly large fraction 
of the excess. These fluxes therefore exclude the contaminating flux 
from the nearby object and are consistent with the spectral energy 
distribution (SED)-derived dust temperature. 

Converting dust temperatures to stellocentric radii is uncertain 
because dust temperature depends on grain size. Typically, the radius 
derived under the assumption of blackbody emission is an under- 
estimate by a factor of up to five”. Thus, we conclude that the dust 
location might, inthe most extreme case, be as far as 1 AU, but not suf- 
ficiently far to explain the 900-day delay between the WISE and optical- 
flux variations, which requires a distance of at least 2 AU. 


Observations 

The beginning of the eclipse was announced? by the ASAS-SN survey, 
which prompted several observing campaigns at optical wavelengths 
andan ALMA observation at Band 7 (programme 2019.A.00040.S). The 
photometric data and passband filters are listed in Extended Data Table 2. 


The All Sky Automated Survey (ASAS)**® is a survey consisting of 
two observing stations, one in Las Campanas, Chile and the other on 
Maui, Hawaii. Each observatory is equipped with two charge-coupled 
device (CCD) cameras using V and I filters and commercial f= 200 mm, 
D=100 mn lenses, although both larger (D = 250 mm) and smaller 
(50-72 mm) lenses were used at earlier times. Most of the data are taken 
with a pixel scale of about 15”. ASAS splits the sky into 709 partially 
overlapping 9° x 9° fields, taking—on average—150 3-min exposures per 
night, leading toa variable cadence of 0-2 frames per night. Depending 
onthe equipment used and the mode of operation, the ASAS limiting 
magnitude varied between 13.5 and 15.5 mag in V and the saturation limit 
was 5.5 to 7.5 mag. Precision is around 0.01-0.02 mag for bright stars 
and below 0.3 mag for the fainter ones. ASAS photometry is calibrated 
against the Tycho Catalogue and its accuracy is limited to 0.05 mag for 
bright, non-blended stars. 

The ASAS-SN’*" consists of five stations around the globe, with 
each station hosting four telescopes with a shared mount. The 
telescopes consist of a 14-cm aperture telephoto lens with a field of 
view of approximately 4.5° x 4.5° and an 8.0” pixel scale. Two of the 
original stations (one in Hawaii and one in Chile) were initially fitted 
with V-band filters, but now these and all the other stations (Texas, 
South Africa and a second station in Chile) observe with g’-band filters 
down to 18 mag. 

The Remote Observatory Atacama Desert (ROAD)* is a fully auto- 
mated telescope located in Chile that obtains nightly photometry in 
Astrodon B, V and I bands for a wide range of astronomical projects. 
It consists of a40-cm //6.8 Optimized Dall-Kirkham and uses a Finger 
Lakes Instruments camera with a 4k x 4k array with pixels of 9 um in 
size. Data are reduced using a custom pipeline and then published on 
the AAVSO website. 

LCOGT is a network of 25 fully robotically operated telescopes 
distributed over seven sites located all around the globe. These 
telescopes are designed to observe transient astronomical events at 
optical and near-infrared wavelengths. LCOGT provides a large variety 
of filter options, but the data we collected are in SDSS g’, r’ andi’ bands. 
All data are automatically processed and calibrated by the BANZAI 
pipeline. The visual companion caused complications in the automatic 
aperture extraction routine of BANZAI; sometimes, correct apertures 
were extracted for both ASASSN-21qj and the nearby star and some- 
times both sources were extracted in one large aperture, often with 
an offset from the true centre of ASASSN-21qj. To correct this, the last 
two stages of the BANZAI routine, aperture extraction and photometry 
calibration, were modified for this specific situation. The calibrated 
magnitudes of all sources in the frames are computed using the default 
BANZAI photometry calibration routine. 

ATLAS is a project that searches for near-Earth asteroids down toa 
magnitude of 19 (ref. 37). Two filters were obtained, the o (orange) and c 
(cyan) filters, respectively. The data consist of two to four photometric 
points observed each night when conditions permitted. Photometry 
with large errors was rejected in a first pass and then the remaining 
observations during a night were averaged and an error based on the 
r.m.s. of these nightly points was calculated. The photometry covers 
the time period in which the collision event occurred. 

The Transiting Exoplanet Survey Satellite (TESS) is a satellite 
designed to survey for transiting exoplanets among the brightest and 
nearest stars over most of the sky. The TESS satellite orbits the Earth 
every 13.7 days ona highly elliptical orbit, scanning a sector of the sky 
spanning 24° x 96° for a total of two orbits, before moving on to the 
next sector. It captures images at 2-s (used for guiding), 20-s (for 1,000 
bright asteroseismology targets), 120-s (for 200,000 stars that are 
probable planet hosts) and 30-min (full-frame image) cadences. The 
instrument consists of four CCD cameras, each with a field of view of 
24° x 24°, with a wide band-pass filter from 600 to 1,000 nm (similar 
to the/.band) and provides high-precision (approximately millimag) 
light curves for stars down to about 14 mag (lo). 
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The NEOWISE mission uses a space-based infrared telescope that has 
been surveying the sky since 2013 at 3.4 and 4.6 um. NEOWISE orbits 
near the Earth’s day-night terminator, scanning rings of the sky near 
about 90° solar elongation and obtains a sequence of observations of 
agiven region of sky every 6 months. The two wavelength channels are 
obtained simultaneously through a beam splitter, allowing for colour 
information to be extracted for each source detected in both bands. 
Detailed descriptions of NEOWISE operations and early results from 
the Reactivation mission” and the standard data processing and data 
characteristics” are available. A single data point for each epoch and 
wavelength is calculated by taking the weighted average of the indi- 
vidual NEOWISE measurements. 

The ALMA data from programme 2019.A.00040.S were downloaded 
and processed through to a measurement set with CASA“. These obser- 
vations were taken on 28 September 2021 (MJD = 59485) and used band 
7, with a mean wavelength of 880 um. No sourceis visible in the default 
archive products and we also detect no source at the expected location 
in CLEAN images. The archive products report anr.m.s. of 17 Wy and 
we measure anr.m.s. of 20 wy ina naturally weighted image. We there- 
fore consider these results as an upper limit of 60 yy. In terms of the 
infrared excess visible in the mid-infrared with WISE, this upper limit 
is not at all constraining but does set limits on emission from cooler 
dust (Fig. 1). 


Light curves 

In this section, we consider the light curves; first, the implications 
of the TESS data for the stellar age and then the implications of the 
ground-based optical light curves for transverse velocity of the occulting 
material and dust grain sizes. 

Stellar rotation and age. TESS data from the Quick-Look Pipeline? ® 
were retrieved from the Mikulski Archive for Space Telescopes (MAST) 
and are shown in Extended Data Fig. 2. The star was observed in sectors 
8,34 and 35. The star was observed soon after the infrared brightening 
and two sequential sectors some time later. 

A periodic signal with a period of approximately 4.3 days is seen, 
with a peak-to-peak amplitude of 2% of the mean flux and the similar 
period can be seen in the two later sectors but they are overwhelmed 
with the first signs of debris from the transiting object. We carry out 
aLomb-Scargle periodogram on S08 and obtain a notable detection 
of P= 4.43 + 0.33 days. A similar period is seen at lower significance in 
the later sectors. Apart from a peak at 1 day, the strongest peak inthe 
ASAS-SN data has a similar period, at 4.1 days. We attribute this modu- 
lation to star spots on the star rotating in and out of our view, and so 
we assert that is the rotational period of the star. Using this rotational 
period, the gyrochronological age*** is calculated to be 300 + 92 Myr. 


Duration of the eclipse and gradient analysis. Figure 2 shows the 
optical light curves. The start of the optical eclipse is seen around MJD 
59350 in the g’ band observations and returns to pre-eclipse levels by 
MJD 59850, giving a total eclipse duration of approximately 500 days. 
The eclipse depth varies as a function of wavelength, whichis discussed 
below. The normalization of the light curves for the ASAS-SN, ATLAS 
and ALLWISE photometry was performed by calculating the 
out-of-transit flux before MJD 58700. The LCOGT and AAVSO photo- 
metry was determined by aligning the photometry of the eclipse with 
the ASAS-SN and ATLAS photometry on a per-band basis. By treating 
the linear changes in flux as owing to the edges of large dust clouds 
crossing the disk of the star, and the largest amount of absorption in 
any of these segments as an estimate for the absorption of the cloud 
(see equations (4.2) and (4.3) in ref. 46), we can determine a robust 
lower limit to the transverse velocity of 7.5 km s” for the material mov- 
ingin front of the star. The method is illustrated in Extended Data Fig. 3 
and the results over the full optical transit are shown in Fig. 2b. 
We convert the magnitude M(t) at time t to a normalized flux 


f(t) =10%©ol25) in which Mo is the mean magnitude outside the 
eclipse. We visually determine turning points in the linearly increasing 
or decreasing photometric flux, fit straight lines to the selected points 
and determine the flux gradient in units of day. A lower bound canbe 
derived for the transverse velocity of the dust, v, by measuring the 
gradient of the light curve and determining what velocity a sharp-edged 
and completely opaque occulter moving across the disk of the star 
would need to make the same gradient. If the dust is on a circular orbit, 
it therefore has to be within 16 AU (equivalent to an orbital period of 
63 years) around the star. The temperature of dust grains at this dis- 
tance would be in the range 100-200 K. 


Dust properties from optical colours. Asis clear from Fig. 2, the pho- 
tometry shows deeper absorption at bluer wavelengths compared with 
redder wavelengths. This wavelength-dependent absorption is typical 
of extinction because of particles with a characteristic size similar to or 
smaller than that of the observed wavelengths. The differences are quan- 
tified in Extended Data Fig. 4, which shows the AAVSO BV/ photometric 
colours asa function of Vmagnitude. For both B- Vand V-/, the colour 
becomes redder as the star dims and the reddening is quantified here 
by the total to selective extinction ratio A,/A, (dashed lines). The values 
are similar to those seen for interstellar extinction”, indicating that 
the dimming is caused by submicrometre-sized dust*®*”’. A further 
constraint from the colours relates to scattering; in protoplanetary disk 
systems that undergo dimming (for example, UX Oritypes), the colour 
initially reddens with dimming but moves back towards the stellar 
colour when the dimming is more than one or two magnitudes (for 
example, ref. 50). This ‘blueing’ is typically interpreted as the relative 
increase in dust-scattered starlight from the disk or envelope as the 
star itself fades”. This behaviour is also seen in Extended Data Fig. 4, 
in which B-V shows marked blueing, whereas V-/ does not. The solid 
lines show an extinction model in which an underlying scattered-light 
component withthe same colour as the star has been added; as the star 
dims, it reddens, but will eventually return to the stellar colour. This 
happens more quickly for B-V because nearly all of the stellar flux in Bis 
blocked and is less pronounced for V-/ because the star is substantially 
less dimmed in /. The fraction of scattered light is relatively large at 
7.5%, implying an important complex of small dust around the star 
by the time the deepest parts of the optical transit occur. This high 
fraction suggests that the impact occurred a substantial fraction of an 
orbit before the optical transit, thus allowing the dust complex time 
to spread around the star. 


Dust-mass estimates 

The SED in Extended Data Fig. 1 gives an estimate of the infrared flux that 
would arise from dust thrown off in the putative collision. We assume 
that the collision occurred on the order of 10 au from the star and hence 
that a dust temperature of 100 K is reasonable. Here 1% of the mass of 
a20 Mean Planet is assumed to be converted entirely into 0.1-um-sized 
grains; this is an optimistic assumption given that the total mass thrown 
offin collisions is on the order 1% or less*” and that this ejected mass is 
constituted of bodies of a range of sizes and all the mass is not present 
as dust. The dust spectrum in Extended Data Fig. 1is approximated as 
a blackbody multiplied by 0.01 x (10/A)'° (with A in um), based on the 
absorption efficiency for silicates”. Changing the grain size tol ym 
yields a similar dust spectrum; although the dust-emitting areais less, 
these grains emit more efficiently. Thus, the mass in small grains detect- 
able with ALMA is about 0.2 Mann. The dust spectrum lies well below 
the WISE measurements and just belowthe ALMA measurement. Even 
in this optimistic case, thermal emission from dust thrown off in the 
collision is therefore not necessarily easily detected. This difference 
is one of surface area; at 600 pc, thousands of square AU of 100-K dust 
is needed for a thermal detection with ALMA, but only a small frac- 
tion of a square AU is needed to notably dim a star (which has radius of 
approximately 0.005 Av). 


Alternative explanations for the observations 

We consider three possible scenarios to explain the observations: 
(1) the initial brightening and later eclipse are two unrelated phenomena; 
(2) the infrared emission and optical transit are both produced by a 
debris disk at about 0.1 AU; (3) we are observing the aftermath of a col- 
lision between super-Earths or mini-Neptunes at a semimajor axis of 
several AU (see Extended Data Fig. 5). We provide here some more detail 
onthe failures of the first two anda preliminary simulation illustrating 
how an impact between two large bodies can produce a large object as 
proposed for our preferred scenario. 

First, the infrared flux increase and the optical dimming might be 
unrelated, coincidental phenomena. For example, one debris disk at 
0.1 AU passively heated to 1,000 K producing the infrared emission and 
another disk further from the star that transited ASASSN-21qj. This 
explanation is unsatisfactory because both infrared flux increases and 
dimming events are rare. Mid-infrared excesses are exceptionally rare 
among main-sequence stars (1:10,000 (ref. 54)) and still uncommon for 
young stars (1:100 (ref. 54)), and no star has previously shown a size- 
able increase starting from no excess. A single case of a disappearing 
mid-infrared excess has been seen, which remains largely unex- 
plained®. Similarly, optical-dimming events are rare for main-sequence 
stars; for example, only one was seen to undergo dust-related optical 
dimming with the Kepler main mission, which observed 150,000 stars 
for 4 years (ref. 56). Both optical and infrared variability are indepen- 
dently less than 1% probabilities for a given star, so for ASASSN-21qj 
to show both by chance is at best a 0.01% probability, and probably 
much lower. 

Another possibility is that both the infrared and optical features are 
produced by a single debris disk at about 0.1 AU from the star. At such 
a close distance, any dust clumps would initially produce periodic 
eclipses on the timescale of days” before being sheared into an 
azimuthally symmetric structure in months, so the occultation of the 
star would need to be related to changes in the vertical structure of a 
near-edge-on post-collision disk, for example, by dynamical ‘stirring’ 
of debris by impact remnants”. The optical depth of the disk must 
then decrease owing to continuing collisional depletion to explain the 
slow return to pre-transit levels of optical flux and gradual decrease in 
infrared flux. Three issues with this model are that: (1) the disk must have 
precisely the right geometry to slowly occult the star as the scale height 
increased and must coincidentally become optically thin; (2) there is 
no apparent change in dust temperature, which would be expected 
as the optical depth decreases and the warmer inner disk becomes 
visible; and (3) notable optical variation is seen 3 years after the putative 
collision, so any initially created clumps would already have sheared 
out. Newer clumps must contribute on the order of 50% of the dust area 
to explain the large variations around MJD 59500 but the infrared flux 
shows a gradual decline, rather than any strong variation that would 
be associated with clump creation. A fourth, but less critical, issue is 
that the inferred clump velocities are not as high as they could be for 
transiting structures at 0.1 AU, which should result in transverse veloci- 
ties of upto100 km s”. 


Smoothed particle hydrodynamics collision simulations 

To provide some more insight into the collision scenario, we performed 
impact simulations using the SWIFT smoothed particle hydrodynam- 
ics (SPH) code” *. Extended Data Fig. 6 shows a collision between two 
25 Mean Planets at 45.77 km s7 (1.4 Ves, escape velocity neglecting the 
atmosphere) at an impact parameter of 0.4 (animpact angle of 23.6°). 
The colliding bodies were 22.5% rock (forsterite), 67.5% water™ and 
10% H/He (ref. 65) by mass. 2.1 x 10° particles were used in the simula- 
tion. To make simulations with high resolution numerically tractable, 
SWIFT imposes a maximum smoothing length that, in effect, imposes 
a minimum density for particles in the simulation (about 30 kg m” for 
the simulation shown here). The bound post-impact material is spread 


over hundreds of Earth radii following the collision, illustrating that 
giant impacts can produce very large post-impact objects. 

However, for the massive and highly extended bodies produced by 
collisions between super-Earths and mini-Neptunes, a large fraction 
of the post-impact body is at the minimum density (green particles in 
lower-right panel in Extended Data Fig. 6). SPH simulations therefore 
probably underestimate the extent of such post-impact bodies and 
further work is needed to fully quantify the size of post-impact bodies 
produced in different impacts. 


Post-impact body-cooling calculations 

How the emission from a post-impact body would evolve with time is 
highly dependent onthe initial mass distribution and thermal state of 
the body, and the balance between radiative cooling, viscous spread- 
ing and mass and angular momentum transport by condensates”??? , 
Given the limitations of SPH simulations (see above), it is not possible 
to accurately determine the initial structure of post-impact bodies in 
the relevant regime. To explore a range of possible evolution pathways, 
we have calculated the evolution of post-impact bodies with different 
power-law surface-density profiles, Z, under the limiting case that 
radiative cooling and condensation of the vapour is the single driver for 
evolution of the structure. We chose a power-law surface-density pro- 
file as it can straightforwardly cover the wide range of surface-density 
profiles expected after super-Earth/mini-Neptune collisions, based 
on those found in impact simulations between lower-mass, terrestrial 
bodies”°°* , The surface-density profile is given by: 


Yo, if nysR, 
Z= ark, if Re< Tey < Remit (1) 
O, if Ky > Remit 


in whichr,, is the distance from the rotation axis; R, is the outer radius 
of aconstant-surface-density central region, roughly analogous to the 
corotating regions seen in Earth-mass synestias””°; 5, is the surface 
density of the central region; £ is the power-law exponent; and Remit iS 
the initial emitting radius. Imposing surface-density continuity gives: 


a= RP. (2) 
We can determine 2, by fixing the mass of the body 


Re Remit 
m-f, Qnty Lodriy J amr, Elri )dri, (3) 


which can be solved to give: 


M, , 
5 3 , if B=-2, 
TIRE + 2R” (InRemit — INR,)] 
Zo= M, ; (4) 
ae , Otherwise. 
nR + Fs Ra RE?) | 


The time taken for a given region of the structure to cool to the point 
that a sufficient fraction of material is condensed for the temperature 
to drop below the condensation buffer (see above)—and thus for the 
emitted flux to drop—is given by 


fyl) 


teooll xy oTt. (5) 


emit 


in whichfis the initial vapour fraction at that radius, lis the latent heat 
of vaporization of the material (here we have taken the limiting case 
of pure water l = 2.256 x 10° (ref. 70), but the addition of silicates could 
make the latent heat much larger), ois the Stefan—-Boltzmann constant 
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and Tomi is the emission temperature. Figure 3 shows example evolu- 
tions of the emission from a post-impact body using this model for 
different parameters (solid lines). Results of a modified model in which 
the initial surface density (by the addition of a constant parameter to 
equation (1)) is forced to be zero at Remix are Shown as dashed lines. 


Data availability 


The datasets generated and analysed during this study are available 
inthe Zenodo repository at https://doi.org/10.5281/zenodo.8344755. 


Code availability 


All the code for the analysis and the generation of all the figures are 
available in a showyourwork” reproducible framework available as a git 
repository at https://github.com/mkenworthy/ASASSN-21qj-collision/. 
The source code and documentation for the SWIFT open-source simu- 
lation code are available from www.swiftsim.com. 
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Extended Data Fig. 1| Spectrum of ASASSN-21qj components. The red 
symbols show the optical and pre-brightening WISE infrared photometry and 
the blue symbols show the post-brightening WISE and ALMA fluxes. Stellar and 
1,000 K components consistent with the pre-brightening and post-brightening 
fluxes are shown. The dashed line shows an estimated cool-component 
spectrum for 0.1-um-sized grains associated with the transiting dust cloud. 
Downward triangles are upper limits. Error bars are shown at lø confidence. 
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Extended Data Fig. 2 | The light curve of ASASSN-21qj from TESS and the 


periodogram of TESS and ASAS-SN photometry. a,b, Photometry of ASASSN- 


21qj from three sectors of TESS. c, Lomb-Scargle analysis of photometry from 
TESS (coloured blue) and from ASAS-SN V-band data from MJD 57420 to 
MJD 58386 (light grey) shows a signal at 4.4 days. At lower frequencies, the 


ground-based photometry shows power and aliasing signals. The TESS signal 
shows an important signal at 4.43 days, anda similar signal is seen inthe 
ground-based data. The longer time baseline in the ASAS-SN data reveals 
substructure in the signal. 
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Extended Data Fig. 3 | Deriving the transverse velocity froma light curve. light-grey vertical lines. b, Gradient of the light curve as a function of time. 
a, ASAS-SN g’ photometry is shown in units of normalized flux. Straight-line c, Transverse velocity derived from the light curve and the gradient of the light 


fits (light-blue lines) are made tothe photometryintheregionsindicatedbythe curve. Error bars are shown at loconfidence. 
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Extended Data Fig. 4 | Blueing of the B-V and V-/ colours during the dimming 


event. Points show AAVSO data and lines show models. a, The Vmagnitude 
versus V-/ colour and b, the Vmagnitude versus V-/ colour. The dashed line is a 
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V-I 


line ofA,/A,for the value shown in the legend and the solid line is a model that 
includes an underlying scattered-light component with s = 7.5% of the stellar 
flux. Error bars are shown at loconfidence. 


Orbital motion 


Dust expands and cools, shearing by 
Keplerian orbital motion and causes 
optical eclipse for 600 days 


t~1000d To Earth 


Remnant forms R~7 solar radii, T~1000K 
photosphere in dust cloud 


t=0 e@ Collision generates dust and hot remnant 
Extended Data Fig. 5 | Sketch of the hypothesis for the observations seen luminosity, generating the 1,000 K infrared emission. Around 1,000 days later, 
towards ASASSN-21qj. At t = 0, the collision occurs, producing a cloud of the expanding cloud crosses the line of sight between the star and the Earth, 


debris that expands and cools. Material close to the remnant is heated by its generating the optical light curve. 
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Bound particles 


Extended Data Fig. 6| Simulations of the formation of a post-impact body. outwards in the initial bodies) and whether they came from the impactor or 
Giant impacts between super-Earths and mini-Neptunescanproducepost-impact target (see top-left panel). The final two panels show just the mass bound to the 
bodies hundreds of Earth radii across, comparable with that requiredtoproduce primary remnant, which has a mass of 48.4 M,,,:,. IN the final panel, particles 
the observed infrared flux. With the exception of the lower-right panel, particles that areat the minimum density imposed by the code are coloured in green. 
are coloured by their material (forsterite, water or aH,-He mixture moving 


Extended Data Table 1| Properties of ASASSN-21qj 


Property Value 
arcrs, bhh mm ss] 08:15:23.301 
ôrcrRs, [dd mm ss] -38:59:23.31 
Ha [mas yr7!] —9.692 + 0.0121 
us [mas yr—+] 7.349 + 0.0121 
w [mas] 1.763 + 0.011! 
RV [km s~+] 25.8 + 3! 
Distance [pc] OTE 
G [mag] 13.371 £ 0.0037 
Gpp [mag] 13.697 + 0.0031 
Grp [mag] 12.882 + 0.0041 
Gpp — Grp [mag] 0.815 + 0.0051 
J [mag] > 12.073 

H [mag] 12.03 + 0.043 
Ks [mag] 11.99 + 0.043 
B (AB) [mag] 14.16 + 0.064 
V (AB) [mag] 13.48 + 0.034 
g' (AB) [mag] 13.77 + 0.034 
r’ (AB) [mag] 13.29 + 0.034 
i' (AB) [mag] 13.20 + 0.084 
Teg [K] 5760 + 10! 
[Fe/H] [dex] —0.23 + 0.014 
logg [logio cms~?] 4.339 + 0.0051 
Rx [Ro] 1.045 
log(Lpor/Lo) [dex] 0.0335 
Mbol [mag] 13.49 0.026 
E(Gpp — Grp) [mag] 0.01 + 0.005° 
Age [Myr] 300 + 927 


'Gaia DR3 (ref. 72), coordinates are J2000 at epoch 2000.0. *Reference 73. °2MASS”, J is an 
upper limit and H and K, are flagged as contaminated. “APASS”, this photometry includes a 
nearby star. °Estimated from an SED fit fixed to Gaia properties. “Estimated using mean stellar 
properties”. "Estimated using rotation period. 
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Extended Data Table 2 | Photometric observations of 
ASASSN-21qj 


Survey name Filter Number of points 


ASASSN V 758 
g’ 3225 
ATLAS c 161 
o 677 
AAVSO B 729 
I 728 
V T22 
LCOGT g 275 
r 224 
i 168 
NEOWISE W1 18 
W2 18 


The number of points are listed per survey and filter. This is the count after initial rejection 
of photometric points with notably large error bars. Further photometric points may have 
been rejected in the different analysis steps. The ROAD observations constitute most of the 
observations from the AAVSO datasets. 


